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ABSTRACT 

An  algorithm  has  been  designed  to  provide  near  real-time  location 
and  classification  information  to  an  operator  of  an  underwater  acoustic 
imaging  system.  The  overall  system  to  detect  objects  buried  up  to 
five  meters  in  unconsolidated  marine  sediments  consists  of  a  spe- 
cialized projector  co-located  with  a  five  meter  receiving  line  array. 
This  system  is  described  in  basic  terms  together  with  the  principles  of 
the  signal  processing  that  extracts  information  from  the  scattered 
acoustic  signals.  Using  simulated  objects,  the  angular  resolution  and 
imaging  performance  of  a  number  of  system  designs  were  evaluated.  To 
obtain  near  real-time  execution,  the  algorithm  was  optimized  and  parti- 
tioned using  parallel,  pipeline,  and  double  buffering  techniques  for 
independent  but  synchronized  operation  on  multiple  microprocessors. 
The  feasibility  of  the  design  approach  was  experimentally  demonstrated 
using  four  single  board  computers  in  a  microcomputer  development  system. 
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I.   INTRODUCTION 

The  specific  situation  which  led  to  this  research  is  the  desire 
to  detect,  locate,  and  classify  objects  buried  in  unconsolidated 
marine  sediments  up  to  a  depth  of  5  meters.  To  accomplish  this  an 
acoustic  imaging  system  is  needed  which  can  be  mounted  on  an  under- 
water vehicle  maneuvering  close  to  the  ocean  floor.  The  problem  is  to 
penetrate  the  soft  sediment  acoustically,  detect  the  scattered  signals, 
and  derive  from  them  information  as  to  the  nature  of  the  scattering 
objects.  With  sufficient  range  and  angular  resolution  it  would  be 
possible  to  adequately  localize  and  classify  an  object.  Furthermore, 
it  is  necessary  to  obtain  this  information  as  rapidly  as  possible: 
ideally  a  display  would  provide  a  "real-time"  representation  of  buried 
objects  to  an  operator  as  a  section  of  sediment  is  being  searched 
acoustically. 

Although  a  number  of  techniques  of  acoustic  imaging  exist  L^l* 
a  system  similar  to  a  high  resolution  sonar  has  been  proposed  which 
would  have  high  resolution  using  very   short  baseband  pulses  and  which 
would  operate  with  dynamic  focussing  deep  in  the  near  field  of  the 
receiving  aperture.  This  system  would  involve  a  specialized  para- 
metric projector  co-located  with  a  receiving  horizontal  line  array  of 
hydrophones.  It  would  be  mounted  on  an  underwater  vehicle  operating  at 
about  5  meters  above  the  ocean  floor,  and  would  provide  acoustic 
penetration  of  the  marine  sediments  to  depths  of  about  5  meters.  The 
objective  of  this  thesis  was  to  design  the  signal  processing  algorithm 
for  such  a  system. 


In  the  course  of  the  development  of  the  algorithm  the  effect  of  a 
number  of  system  parameters  on  angular  resolution  and  imaging  perform- 
ance was  evaluated  for  a  number  of  possible  implementations  of  the 
acoustic  system.  This  led  to  major  decisions  as  to  the  form  of  the 
algorithm  and  the  possible  acoustic  imaging  system. 

Since  an  acoustic  imaging  system  mounted  on  a  maneuvering  under- 
water vehicle  would  be  required  to  provide  an  operator  with  a  near 
real-time  display  of  buried  objects,  the  algorithm  was  redesigned  for 
increased  execution  efficiency  and  implementation  on  multiple  micro- 
processors. Optimization  of  the  algorithm  by  removing  time  consuming 
operations,  together  with  incorporation  of  parallel  and  pipeline 
processes  and  double  buffering  concepts  were  the  approaches  taken. 
Experiments  were  conducted  using  a  number  of  single  board  computers  to 
test  and  demonstrate  the  feasibility  of  near  real-time  execution  of  the 
algorithm. 

Section  II  of  the  thesis  describes  the  acoustic  imaging  system  and 
the  basis  for  the  signal  processing  algorithm.  The  details  of  the 
algorithm  and  the  angular  resolution  and  imaging  performance  evaluations 
are  explained  in  Section  III.  The  design  of  the  algorithm  for  near 
real-time  execution  on  multiple  microprocessors  and  the  experimental 
demonstration  of  feasibility  are  contained  in  Section  IV.  The  resulting 
conclusions  are  then  given  in  Section  V,  followed  in  Section  VI  by 
recommendations  for  future  work  that  could  produce  a  special  purpose 
processor  for  the  acoustic  imaging  system. 
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II.   THE  PROPOSED  ACOUSTIC  IMAGING  SYSTEM 

A.  INTRODUCTION 

An  acoustical  imaging  system  mounted  on  an  underwater  vehicle  will 
consist  essentially  of  three  major  parts:  a  sophisticated  projector,  a 
receiving  array,  and  a  signal  processor  and  display  system.  The  thrust 
of  this  thesis  has  been  to  develop  the  signal  processing  algorithm.  To 
understand  the  underlying  concepts,  the  potential  form  of  the  projector 
and  receiving  array  is  briefly  described,  followed  by  the  proposed 
algorithm. 

B.  TRANSMISSION 

Acoustic  waves  can  penetrate  sediment  to  useful  depths  with 
acceptable  absorption  loss  if  the  frequency  is  low  enough;  less  than 
about  20  kHz.  Very  short  acoustic  pulses  (50  to  100  microseconds)  are 
required  for  adequate  range  and  angular  resolution.  For  the  concept 
that  has  been  developed  in  this  thesis  it  has  become  apparent  that  it 
would  be  necessary  to  have  effectively  a  point  source  for  the  projector. 
Baseband  pulse  generation  using  a  parametric  source  [2]  appears  to  be  a 
promising  technique  for  this  application.  Such  a  projector  would  form 
a  pencil  beam  in  the  collimated  region  of  its  near  field  that  could  be 
scanned  to  produce  a  tomographic  image  of  a  region  in  the  sediment. 

C.  RECEPTION 

1 .  Physical  Configuration 

The  scattered  acoustical  energy  would  be  received  by  a  horizon- 
tal line  array  (mounted  on  the  underwater  vehicle)  of  a  large  number  of 
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equally  spaced  hydrophones.  The  projector  would  be  located  at  the  center 
of  the  line  array.  The  length  of  the  array  would  be  a  compromise  between 
the  need  for  a  large  aperture  (on  the  order  of  10  meters  for  angular 
resolution)  and  the  physical  reality  of  deploying  and  maneuvering  an 
underwater  vehicle  with  the  array  attached  (5  meters  probably  the 
maximum). 

Moving  over  the  area  of  interest  the  underwater  vehicle  would 
maintain  an  altitude  of  approximately  5  meters  above  the  sediment.  As 
the  desired  depth  of  penetration  into  the  sediment  is  also  5  meters  the 
total  range  of  interest  lies  between  5  and  10  meters  from  the  center  of 
the  array.  A  60  degree  sector  (between  array  bearings  of  -30  and  +30 
degrees  from  the  normal  to  the  line  array)  provides  horizontal 
coverage  of  10  meters  in  the  sediment  at  10  meters  range  from  the 
center  of  the  line  array  (Figure  1). 
2.  Wavefront  Curvature 

For  ranges  between  5  and  10  meters,  operation  of  the  acoustical 
imaging  system  is  deep  in  the  near  field  of  the  line  array.  The  far 
field  begins  at  a  range  of  D^/A  where  D  is  the  diameter  of  the  aperture 
and  X    is  the  wavelength  [Ij  .   Using  a  50  microsecond  pulse  and  a  5 
meter  line  array  this  is  on  the  order  of  300  meters.  The  near  field 
operation  is  evident  by  the  marked  curvature  of  the  wavefronts  of  the 
acoustic  energy  as  they  arrive  at  the  line  array  after  scattering  from 
points  within  the  area  of  interest.  Since  the  curvature  depends  on  the 
range  to  a  scattering  point  an  imaging  processor  requires  some  form  of 
dynamic  focussing. 
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Figure  1.     The  physical  relationship  between  the 
five  meter  line  array  of  hydrophones 
and  the  sector  of  interest. 
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Figure  2  depicts  four  such  wavefronts  (that  originated  at  the 
indicated  bearings  and  ranges)  as  a  function  of  time  of  arrival  and 
position  along  a  5  meter  line  array.  Using  1500  meters  per  second  as 
the  soundspeed,  the  time  of  the  earliest  arrival  of  interest  is  6.220 
milliseconds  after  pulse  transmission,  and  the  latest  is  14.304  milli- 
seconds. 

3.  The  "Trace"  of  a  Point  Scatterer 

The  line  array  consists  of  a  large  number  of  equally  spaced 
discrete  hydrophones  (each  considered  as  a  point)  which  detect  the 
pressure  amplitude  of  an  arriving  wavefront.  The  time  sampling  of  the 
outputs  of  16  hydrophones  in  a  5  meter  line  array  resulted  in  a  quanti- 
zation of  the  wavefront  arrival  times  as  shown  in  Figure  3  (in  which 
an  acoustic  pulse  width  of  50  microseconds  and  sampling  ewery   50  micro- 
seconds, beginning  6  milliseconds  after  pulse  transmission,  was  used). 
For  example,  considering  the  wavefront  originating  at  a  scattering 
point  located  at  a  bearing  of  zero  degrees  and  a  range  of  five  meters, 
the  hydrophone  sampling  produces  the  following  time  indices  for  the 
wavefront  arrivals:  21,  19,  17,  16,  15,  14,  13,  13,  13,  13,  14,  15, 
16,  17,  19.  21. 

The  sequence  of  time  Indices  representing  the  arrival  of  a 
wavefront,  n_(r,9)  *  (n,  ,n2...n.)  for  K  sensors,  is  unique  for  the 
scattering  element  at  a  range  r  and  bearing  9:  this  sequence  is  termed 
the  "trace"  of  that  particular  point  (related  to  the  concept  of  the 
"Time  Delay  Trace  Function"  [3]  ).  The  traces  for  all  points  at  the 
range  and  bearing  increments  in  the  area  of  interest  can  be  pre-calcula- 
ted  and  stored  for  later  retrieval . 
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Figure  2.     The  arrival   times  of  four  wavefronts  as 

a  function  of  position  along  a  five  meter 
array. 
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Figures.     Quantization  of  wavefront  arrival   times. 
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D.   SIGNAL  PROCESSING 

1 ,  The  Basic  Concept 

In  developing  the  algorithm  two  situations  were  examined: 

a.  Insonification  of  the  entire  sector  by  a  short  pulse 
(broad  beam  in  the  plane  of  the  line  array  but  very  narrow  beam  in  the 
plane  perpendicular  to  the  line  array)  followed  by  processing  returns 
from  the  entire  sector,  and 

b.  Transmission  of  the  acoustic  pulse  in  a  narrow  beam  and 
sweeping  across  the  sector  with  signal  processing  occurring  at  each 
bearing  increment. 

Essentially,  after  a  suitable  gate  time  after  pulse  trans- 
mission, each  sampled  hydrophone  output  provides  a  time  series  record 
of  the  scattered  acoustic  amplitudes  arriving  at  the  line  array. 
Storing  the  results  (after  A/D  conversion)  at  each  sample  time,  n,  of 
all  K   hydrophone  outputs  causes  the  memory  to  contain  a  time  represen- 
tation of  the  scattered  acoustic  pressure  field,  D(n,k),  k=l..K,  n=l .  .N 
where  N  is  the  maximum  sample  time.  The  presence  of  a  scattering  point 
will  be  indicated  by  the  presence  in  memory  of  signal  amplitudes 
occurring  at  the  appropriate  times  in  each  hydrophone  record  (a  pattern 
representing  the  wavefront  curvature). 

The  sequence  of  time  indices  which  uniquely  represent  a  scat- 
tering point  are  available  as  the  pre-calculated  trace  for  that  point. 
An  estimator  of  the  presence  of  a  pattern  is  the  summation  of  the  signal 
amplitudes  at  the  appropriate  time  indices  or  addresses  in  each  of  the 
hydrophone  records.  For  example,  if  the  signal  amplitude  at  the 
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wavefront  (from  a  point  at  r,9)  was  considered  unity,  then  the  total 
summation  along  ]i(r,9)  would  equal  K,  the  total  number  of  sensors: 

K 
A(r,9)  =X  D(n^.k)  =  K 
k=l 
where  n^^  are  elements  of  n^(r,9). 

As  a  result  of  noise  and  errors  in  geometry  the  sum  would  be  less  than 
K  but  would  nevertheless  have  a  maximum  if  a  scattering  point  was 
present. 

In  summary,  what  is  required  is  to  form  in  memory  a  represen- 
tation of  the  scattered  acoustic  field,  search  through  it  systematically 
for  wavefront  patterns  by  summation  along  pre-calculated  traces  for  the 
ranges  and  bearings  of  interest,  and  recording  the  resultant  total 
amplitudes  as  a  function  of  range  and  bearing.  A  decision  process 
would  then  follow  to  determine  which  amplitudes  represented  point 
scatterers,  as  opposed  to  solely  noise,  and  then  these  could  be  dis- 
played forming  the  acoustical  image  of  the  insonified  region  in  terms 
of  range  and  bearing  or  rectangular  coordinates. 
2.  Data  Structures 

To  visualize  the  signal  processing  operation  it  is  useful  to 
discuss  it  in  terms  of  the  data  structures  involved:  the  "Indata" 
array,  the  "Trace"  matrix,  and  the  "Amplitude"  array  as  shown  in 
Figure  4. 

The  Indata  array  contains  the  time  record  of  the  scattered 
acoustic  field.  It  can  be  visualized  as  two  dimensional  in  memory: 
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Figure  4.   The  basic  data  structures 
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i.e.  D(n,k)  where  n  is  the  time  index,  and  k  is  the  hydrophone  number. 
Each  row  would  be  formed  by  the  simultaneous  sampling  of  all  hydrophones 
at  intervals  equal  to  the  pulse  width  of  50  microseconds,  beginning 
after  a  gate  time  of  6  milliseconds,  and  storing  the  results.  This 
would  continue  for  167  samples  until  the  last  possible  arrival  time  of 
a  signal  of  interest  (time  indices  0  to  166). 

The  Trace  matrix,  T(9,r,k)  contains  the  traces  n_(r,9)  for 
points  in  th^  area  of  interest  as  identified  by  the  intersections  of  all 
range  and  bearing  increments.  It  can  be  visualized  as  three  dimen- 
sional where  the  first  dimension  is  that  of  bearing,  the  second  is  that 
of  range,  and  the  third  is  the  hydrophone  number.  For  a  particular 
range  and  bearing,  the  elements  in  the  third  dimension  contain  the  time 
index  sequence  n^(r,9)  that  is  the  pre-calculated  trace:  i.e.  nj^  = 
T(9,r,k). 

The  dimensions  for  the  Amplitude  array  are  range  and  bearing. 
Each  element,  A(r,9),  is  the  total  sum  of  the  contents  of  those 
elements  from  the  Indata  array  identified  by  the  trace  n^(r,9): 

K 
i.e.  A(r,9)  =^  D(n,^,k) 

k=l 

where  n.  are  elements  of  ji(r,9) 

Scattering  points  in  the  acoustically  scanned  region  will  be  evident 
by  larger  amplitude  elements  in  A(r,9)  and  the  distribution  of  these 
will  comprise  an  image  of  any  objects  present. 
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III.   THE  DEVELOPMENTAL  ALGORITHM 

A.  INTRODUCTION 

The  signal  processing  algorithm  was  initially  developed  using  a 
Hewlett-Packard  9845B  desktop  minicomputer  (Figure  5).  This  machine 
uses  HP  enhanced  BASIC  [4]  and  was  equipped  with  192K  of  storage,  two 
peripheral  flexible  disk  drives  (HP9885M  and  HP9885S),  a  Graphics  ROM 
(98437B),  and  a  CRT  Graphics  Memory  (98470B).  This  proved  to  be  a  most 
convenient  tool  for  the  development  of  the  algorithm,  evaluation  of 
the  influence  of  the  major  parameters  involved,  and  for  the  graphical 
display  of  simulated  objects,  all  of  which  led  to  major  conclusions  and 
decisions  as  to  the  form  of  the  algorithm  and  the  possible  physical 
acoustic  imaging  system. 

B.  DESCRIPTION  OF  THE  ALGORITHM 
1 .  Parameters 

The  parameters  involved  in  the  algorithm  (Appendix  A)  were 
those  which  describe  the  physical  implementation  of  an  underwater 
acoustic  imaging  system  and  the  area  or  sector  of  interest  within  the 
sediment  that  is  to  be  insonified.  These  parameters  were  varied  to 
evaluate  their  influence  on  the  performance  of  the  acoustical  imaging 
system. 

The  lengths  of  the  horizontal  line  array  used  in  the  develop- 
mental algorithm  were  10  meters  and  5  meters;  indicative  of  the  desire 
for  the  largest  aperture  possible  for  angular  resolution,  and  the  most 
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probable  realistic  maximum  for  a  physical  structure  to  be  attached  to 
a  maneuvering  underwater  vehicle.  Equally  spaced  hydrophones  along 
the  line  array  numbering  16  and  32  were  used,  as  were  pulse  widths  of 
100  and  50  microseconds.  The  sound  speed  was  taken  as  1500  meters 
per  second  throughout  the  transmission  path  of  the  acoustic  pulse. 

The  sector  of  interest  within  the  sediment  was  defined  as 
lying  between  the  ranges  of  5  and  10  meters  from  the  center  of  the  line 
array  and  between  bearings  of  -30  and  +30  degrees  from  the  normal  to  the 
line  array.  Range  increments  of  0.1  meters,  and  bearing  increments  of 
two  degrees  were  chosen.  This  corresponded  to  51  range  indices  from  0 
to  50,  and  31  bearing  indices  from  -15  to  15.  Bearing  increments  of 
one  degree  were  used  in  the  evaluation  of  array  bearing  resolution. 
2.  Calculation  of  Traces 

The  trace,  _n(r,9),  for  a  scattering  point  at  r,6  consists  of 
the  time  indices  representing  the  arrival  of  the  scattered  wavefront 
at  each  of  the  sensors  in  the  line  array.  With  the  projector  at  the 
center  of  the  line  array,  the  distance  from  the  projector  to  each 
sensor  was  calculated.  A  simple  trigonometric  relation  was  used  to 
calculate  the  travel  time  for  an  acoustic  pulse  from  the  projector  to 
a  point  P(r,9)  and  then  to  the  k'th  sensor. 

The  geometry  of  two  cases  is  illustrated  in  Figure  6  and 

Figure  7.  The  first  case.  Figure  6,  applied  when  the  k'th  sensor  was 

past  the  array  midpoint;  i.e.  k  >  K/2  where  K  is  the  number  of  sensors 

in  the  line  array.  For  this  situation  the  trigonometrical  relation 

providing  the  distance  d|^  from  the  k'th  sensor  to  the  point  P(r,9)  in 

terms  of  the  range  r,  bearing  9,  and  sensor  to  projector  distance  bj^,  is 

as  f ol 1 ows  : 
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line  array 


P(r,9) 


Figure  6.       The  geometry  for  k  >  K/2 


line  array 


P(r.8) 


Figure  7.   The  geometry  for  k  ^  K/2. 
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dj^2  =  b|^2  +  y.2  .  2  b|^  r  cos(90-e) 

Similarly,  for  the  case  k  ^  K/2  (Figure  7),  the  distance  dj^  is 
obtained  from 

d|^2  =  b,^2  +  p2  -  2  b,^  r  cos(90+9) 

The  total  travel  distance  for  a  scattered  pulse,  r  +  d]^,  divided  by 
the  sound  speed,  gives  the  time  that  the  leading  edge  of  the  pulse 
arrives  at  the  k'th  sensor  after  transmission.  The  period  of  6  milli- 
seconds is  subtracted  from  the  travel  time,  which  is  equivalent  to 
starting  data  acquisition  after  gating  out  the  first  4.5  meters  in 
range  from  the  center  of  the  line  array. 

The  integer  portion  of  the  quotient  of  this  time  and  the 
sampling  interval  (equivalent  to  the  pulse  width)  provides  the  time 
index  n;^.  Computed  for  each  of  the  K  sensors  this  results  in  the 
sequence  n^(r,9)  =  {n^^  where  k  =  1..K.  This  computation  was  repeated 
at  all  ranges  for  each  of  the  positive  bearings.  Symmetry  about  the 
normal  to  the  line  array  permits  the  calculations  of  travel  times  for 
positive  bearings  to  suffice.  Negative  bearings  use  the  mirror  image 
of  the  trace  sequence;  i.e.  n^(r,-9)  =  {n^  where  k  =  K..1. 
3.  Simulation  of  Objects 

In  order  to  develop  the  algorithm  and  in  the  absence  of  real 
hydrophone  inputs,  sampling,  and  I/O  operations  it  was  necessary  to 
simulate  the  Indata  array  D(n,k).  An  idealized  situation  consisting  of 
zero  noise  and  unity  signal  amplitude  at  the  wavefront  of  a  scattered 
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pulse  was  considered.  In  other  words  no  attempt  was  made  to  model  a 
noise  environment  or  the  reflectivity  of  a  buried  object. 

A  point  scatterer  was  simulated  by  using  the  pre-calculated 
trace,  as  stored  in  the  Trace  matrix  T(e,r,k),  to  identify  the  time 
indices  n  of  D(n,k)  that  represented  the  wavefront  arrival  at  each 
sensor.  Unity  signal  amplitude  was  inserted  into  these  elements.  For 
the  narrow  acoustic  beam  case  only  a  single  point  scatterer  was  assumed 
to  exist  at  each  range  increment  along  the  bearing  of  the  beam.  It  was' 
more  complicated  in  the  case  where  the  entire  sector  was  assumed  insoni- 
fied.  In  this  case  the  acoustic  field  would  be  composed  of  the  inter- 
ference pattern  of  numerous  wavefronts  if  the  sector  contained  a 
distributed  object  composed  of  a  number  of  point  scatterers. 

For  the  latter  case  two  input  situations  from  the  sensors  were 
simulated:  a  nonlinear  input,  and  a  linear  input.  The  nonlinear 
situation  would  arise  if  the  outputs  from  the  hydrophones  were  ampli- 
fied, passed  through  a  thresholding  comparator  circuit,  and  a  binary 
input  presented  to  the  computer  I/O.  The  unity  amplitude  mentioned 
above  simulated  this  situation.  If  the  sampled  outputs  from  the 
sensors  were  subjected  to  A/D  conversion  prior  to  storage  in  memory, 
then  a  more  accurate  or  linear  representation  of  the  actual  pressure 
field  would  be  available.  This  was  simulated  by  the  summation  of  the 
wavefront  amplitudes  in  the  appropriate  elements  of  D(n,k)  as  each  of 
the  points  comprising  a  distributed  object  was  considered. 

The  angular  resolution  of  the  imaging  system  was  evaluated  by 
simulating  discrete  points  in  the  area  of  interest.  Imaging  perform- 
ance was  studied  using  simulated  objects  composed  of  a  number  of 
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scattering  points.  The  resulting  Indata  arrays  were  then  operated 
upon  by  the  algorithm  to  form  the  Amplitude  array  which  represented  the 
acoustic  image  of  the  insonified  sector. 
4.   Forming  the  Amplitude  Array 

The  magnitude  of  each  element  in  the  Amplitude  array,  A(r,9), 
is  an  estimator  of  the  presence  of  a  point  scatterer  at  the  range  r, 
and  the  bearing  9.  Each  element  results  from  the  sum  of  those  elements 
in  the  Indata  array,  D(n,k),  identified  by  the  trace  r[(r,9)  stored  in 
the  Trace  matrix,  T(9,r,k). 

For  the  case  of  the  entire  sector  being  insonified,  the  Ampli- 
tude array  was  formed  by  summing  along  the  traces  for  all  ranges  at 
each  bearing  in  turn  for  the  entire  sector.  For  the  narrow  beam  case, 
only  the  traces  corresponding  to  all  ranges  at  the  beam  bearing  were 
utilized,  followed  by  the  acquisition  of  a  new  Indata  array  for  the  next 
bearing.  Then  the  traces  for  all  ranges  at  that  bearing  were  used  to 
guide  the  summation  process,  continuing  until  the  narrow  beam  had 
scanned  the  entire  sector  and  the  contents  of  the  Amplitude  array 
represented  the  presence  of  any  point  scatterers  within  the  sector. 
5.  Graphics  Display 

To  evaluate  angular  resolution  and  imaging  performance  of  the 
proposed  system  the  graphics  capability  of  the  HP9845B  was  utilized. 
It  was  possible  to  draw  graphs  representing  the  received  amplitudes 
from  point  scatterers  and  thereby  obtain  the  system's  angular  resolu- 
tion. Both  a  polar  or  range-bearing  grid  and  a  rectangular  grid  were 
devised  to  provide  a  means  of  referencing  a  displayed  image.  It  was 
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decided  that  the  rectangular  grid  would  be  more  practical  and 
meaningful  to  an  operator  of  an  acoustic  imaging  system. 

Due  to  the  choice  of  0.1  meter  increments  in  range  and  two 
degree  increments  in  bearing,  a  single  point  scatterer  was  represented 
by  graphically  "filling-in"  an  area  on  the  screen  equivalent  to  these 
dimensions  and  centered  at  the  range  and  bearing  of  the  point.  This 
was  termed  a  "target"  in  the  developmental  algorithm.  A  threshold 
was  used  to  determine  when  the  magnitude  of  an  element  in  the  Ampli- 
tude array  indicated  a  target  present.  An  object  was  simulated  by 
using  a  number  of  these  discrete  points  which  the  graphics  display 
then  presented  as  a  number  of  adjacent  targets.  This  gave  the  image 
outline  a  stepped  appearance  rather  than  the  smooth  lines  of  a  real 
object. 

C.   RESULTS 

1 .  Range  Resolution 

The  use  of  an  extremely  short  pulse  of  50  microseconds  has 
inherent  range  resolution  much  greater  than  required  or  utilized  in  the 
development  of  the  algorithm.  It  was  decided  that  a  range  resolution 
of  0.1  meters  was  adequate  to  localize  and  classify  large  buried 
objects.  Accordingly,  it  was  the  range  increment  chosen  for  the  pre- 
calculation of  all  the  traces,  and  the  subsequent  processing  and 
display  of  the  acoustic  image.  An  additional  consideration  was  the 
memory  size  available  to  store  the  Trace  matrix. 
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2.  Bearing  Resolution 

For  the  case  where  the  entire  sector  is  insonified  by  one  pulse, 
the  effects  of  a  number  of  parameter  changes  on  the  angular  resolution 
of  the  acoustic  imaging  system  were  obtained.  The  bearing  resolution 
was  determined  by  the  ability  of  the  system  to  distinguish  between  two 
discrete  points  at  the  same  range  but  at  slightly  different  bearings. 
The  two  points  were  simulated  as  being  in  the  vicinity  of  a  bearing  of 
zero  degrees,  and  then  a  bearing  of  thirty  degrees,  and  at  ranges  of 
five  meters  and  ten  meters. 

The  algorithm  processed  the  resulting  Indata  array  in  each 
case  and  the  bearing  resolution  was  determined  from  the  resulting 
amplitude  distribution  in  bearing  at  the  range  of  the  points.  Figure 
8  gives  an  example  of  the  amplitude  distribution  for  two  points  at  a 
number  of  bearing  separation  angles.  The  amplitude  distributions  were 
superimposed  as  the  bearing  separation  between  the  points  was  increased 
from  three  degrees  to  six  degrees  in  one  degree  increments.  The  separa- 
tion in  bearing  that  resulted  in  3db  and  6db  dips  between  the  peaks  in 
the  amplitude  distribution  were  recorded  and  are  shown  in  Table  1. 

This  was  done  for  both  the  nonlinear  input  case  and  the  linear 
input  case,  and  for  sixteen  and  thirty-two  hydrophones  in  the  line 
array.  It  was  observed  that  little  difference  in  bearing  resolution 
resulted  as  could  be  expected  for  the  idealized  zero  noise  situation 
that  was  being  simulated. 
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3.  Simulated  Images 

Three  large  objects,  each  considered  as  being  composed  of  a  num- 
ber of  point  scatterers,  were  simulated  in  the  sector  or  the  field  of 
view  of  the  system.  One  was  small  and  rounded  as  if  it  were  a  cross 
section  of  a  long  cylindrical  object.  The  other  two  were  long  narrow 
objects;  one  was  in  a  horizontal  position,  and  the  other  at  an  inclined 
aspect  angle.  Figure  9  shows  the  positions  and  orientations  of  the 
simulated  objects. 

For  the  case  of  a  pulse  insonifying  the  entire  sector  and 
with  nonlinear  inputs  from  the  receiving  hydrophones,  the  maximum 
amplitude  for  any  element  in  the  Amplitude  array  was  equal  to  the  number 
of  hydrophones.  The  threshold  for  the  display  of  the  resulting  targets 
was  varied  to  see  the  spatial  distribution  (in  range  and  bearing)  of 
the  amplitudes  of  the  image  elements.  Figure  10  shows  the  smearing  of 
the  objects  resulting  from  the  lack  of  angular  resolutions;  i.e.  the 
skirts  of  the  amplitude  distribution  in  bearing  of  each  scattering 
point  were  being  detected. 

Similarly  Figure  11  shows  the  results  of  having  linear  inputs 
from  the  hydrophones.  The  largest  amplitude  elements  were  considerably 
larger  than  the  nonlinear  case,  and  the  contours  of  the  amplitude  dis- 
tribution in  the  image  were  much  steeper.  It  was  observed  that  the 
size,  shape  and  aspect  angle  of  the  objects  had  a  substantial  effect 
on  the  resulting  amplitude  distribution  in  the  image.  This  would  be 
beneficial  for  the  identification  or  classification  function  of  an 
acoustic  imaging  system  particularly  if  it  provided  an  operator  con- 
trolled threshold  for  the  display  or  if  the  intensity  of  the  display 
depicted  the  amplitudes  of  the  image  elements. 
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For  the  case  of  a  narrow  two  degree  beam  width  system  that 
transmits  a  pulse  at  each  bearing  increment,  the  resulting  image  after 
the  sector  was  completely  scanned  is  shown  in  Figure  12.  It  was  ob- 
served that  (as  expected)  the  image  exactly  duplicated  the  simulated 
objects. 

D.   CONCLUSIONS 

The  developmental  algorithm  was  used  to  study  some  of  the  effects 
of  the  system  parameters  on  bearing  resolution,  and  the  imaging  per- 
formance of  possible  implementations  of  an  acoustic  system  that  used 
the  proposed  signal  processing  algorithm.  It  was  concluded  that  the 
acoustic  imaging  system  should  be  a  narrow  beam  system  that  scans  the 
sector  by  transmitting  a  very  short  pulse  at  each  bearing  increment, 
provides  linear  inputs  to  memory  from  the  hydrophone  line  array,  and 
provides  an  operator  controlled  threshold  for  the  image  display.  Al- 
though the  benefits  of  the  latter  two  were  not  evident  in  the  simulated 
idealized  narrow  beam  images,  it  is  felt  that  they  would  be  significant 
in  a  real  system  where  noise  and  objects  with  different  scattering  or 
reflective  properties  were  present.  This  is  also  the  basis  for 
requiring  a  line  array  as  a  large  receiving  aperture,  together  with  the 
proposed  signal  processing  algorithm,  since  objects  with  specular 
reflective  properties  would  not  necessarily  scatter  the  acoustic  energy 
in  the  direction  of  the  projector. 
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IV.        THE  NEAR  REAL-TIME  ALGORITHM 

A.  INTRODUCTION 

An  acoustic  imaging  system  mounted  on  a  maneuvering  underwater 
vehicle  would  be  required  to  provide  an  operator  with  a  near  real- 
time display  of  any  objects  buried  in  the  sediment  as  the  vehicle 
conducts  a  search  pattern  over  the  ocean  floor.  Although  the  signal 
processing  algorithm  that  was  developed  was  basically  simple  computa- 
tionally, it  was  nevertheless  very  time  consuming  and  would  have  been 
inadequate  in  an  operational  environment.  Consequently,  it  was  re- 
designed for  increased  execution  efficiency  and  implementation  on 
multiple  microprocessors  while  maintaining  the  basic  signal  processing 
approach  that  has  been  described.  Optimization  of  the  algorithm  by 
removing  time  consuming  operations,  together  with  incorporation  of 
parallel  and  pipeline  processes  and  double  buffering  concepts  were  the 
approaches  taken.  Experiments  were  conducted  using  a  number  of  single 
board  computers  to  test  and  demonstrate  the  feasibility  of  the  rede- 
signed algorithm  for  near  real-time  execution. 

B.  THE  TIME  PROBLEM 

An  underwater  vehicle  acoustically  scanning  the  sediment  and  moving 
at  1  meter  per  second  would  probably  require  a  complete  sector  scan  at 
least  once  per  second  as  a  minimum  adequate  display  rate;  i.e.  a  two 
degree  beam  width,  approximately  0.34  meters  wide  at  10  meters  range, 
would  scan  a  60  degree  sector  perpendicular  to  the  direction  of  travel 
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once  every  meter  of  forward  travel.  At  each  bearing  increment  data 
acquisition  would  begin  6  milliseconds  after  pulse  transmission  and 
continue  for  8.5  milliseconds.  Thus  the  processing  of  the  Indata 
array  can  begin  14.5  milliseconds  after  pulse  transmission,  and  would 
need  to  be  completed  in  about  15  milliseconds  for  each  of  the  bearing 
increments  if  a  scan  were  to  be  completed  in  about  one  second.  It  was 
obvious  that  the  developmental  algorithm  would  require  considerable 
changes  for  an  efficient  microprocessor  implementation. 

C.   ALGORITHM  OPTIMIZATION 

The  algorithm  was  first  examined  with  the  objective  of  removing 
time  consuming  operations  such  as  multiplication.  These  occurred  in  the 
addressing  of  the  multi -dimensional  arrays  in  the  algorithm  which  were 
stored  in  row  major  order.  Four  additional  one  dimensional  arrays  were 
introduced  which  provided  a  table  lookup  of  the  row  addressing.  For 
example,  the  addresses  for  the  bearing  dimension  of  the  Trace  matrix, 
T(Q,r,k),  were  stored  in  B_array(b),  and  those  for  the  range  dimension 
in  R_array(r).  Instead  of  the  addressing  for  a  single  element  of 
T(9,r,k)  requiring  three  additions  and  two  multiplications,  it  was 
addressed  using  T{B_array(b)  +  R_array(r)  +  k);  i.e.  five  addition 
operations.  This  was  done  for  each  of  the  multi-dimensional  arrays,  the 
Indata  array  D(n,k),  the  Trace  matrix  T(9,r,k),  and  the  Amplitude  array 
A(r,9),  thus  reducing  the  data  structures  in  the  algorithm  to  strictly 
one  dimensional  arrays. 

With  the  intention  of  testing  the  feasibility  of  a  microprocessor 
implementation  of  the  algorithm  using  a  number  of  Intel  SBC  80/20  Single 
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Board  Computers  (Intel  8080  CPU)  [s]  ,  the  algorithm  (Appendix  B)  was 
written  using  the  high  level  language  PL/ 1-80  [6^  .  Executing  the 
program  on  one  SBC  80/20  it  was  found  that  after  the  pre-calculation  of 
the  traces  in  T(9,r,k)  the  signal  processing  aspect  of  the  algorithm 
required  6.45  seconds  to  complete  one  sector  scan  (approximately  210 
milliseconds  at  each  bearing,  and  not  including  the  time  required  for 
pulse  propagation  and  data  acquisition). 

D.   ALGORITHM  PARTITIONING 

The  structure  of  the  signal  processing  algorithm  was  examined  in 
order  to  identify  functions  that  could  be  accomplished  by  parallel  and 
pipelined  processes  executing  on  a  number  of  microprocessors  for  in- 
creased time  efficiency.  The  concept  of  double  buffering  was  also 
incorporated  to  reduce  execution  time. 

1 .  Parallel  Processing 

The  processing  carried  out  at  each  of  the  range  increments  on 
the  Indata  array,  for  the  bearing  of  the  transmitted  pulse,  is  indepen- 
dent of  any  of  the  other  ranges.  This  is  the  summation  that  results 
in  an  elanent  of  the  Amplitude  array  A(r,8).  It  was  decided  to  parti- 
tion this  function  into  separate  processes  that  were  identical  except 
for  the  scope  of  the  range  indices.  These  were  called  the  Sum_ampli- 
tude  processes  and  were  executed  in  parallel  on  a  number  of  micropro- 
cessors. Using  three  microprocessors  for  example,  process  S,  performed 
summations  along  the  traces  for  the  range  indices  0  to  16,  process  $« 
used  range  indices  17  to  33,  and  process  S^  covered  range  indices  34  to 
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50.  In  this  manner  all  ranges  were  processed  In  one  third  the  time 
that  was  required  by  a  single  process  executing  all  range  increments. 

2.  Pipeline  Processing 

The  function  of  data  acquisition  is  called  the  Indata  process 
and  was  seen  to  be  amenable  to  pipeline  processing  with  the  Sum_ampli- 
tude  processes.  This  process  could  control  the  timing  of  pulse  trans- 
mission, sampling  the  hydrophone  outputs,  and  storing  the  results  into 
two  buffer  arrays.  The  Indata  array  of  the  developmental  algorithm 
became  two  arrays  called  the  A_indata  array  and  the  B_indata  array. 
After  the  A_indata  array  was  filled  with  the  data  from  a  particular 
bearing,  the  Sum_amplitude  processes  were  permitted  to  begin  their 
operations  on  the  A_indata  array.  Simultaneously,  the  Indata  process 
proceeded  to  fill  the  B_indata  array  with  the  data  from  the  next 
bearing  increment.  When  the  Sum_amplitude  processes  had  finished  with 
the  A_indata  array,  they  processed  the  B_indata  array.  The  Indata 
process  could  then  utilize  the  A_indata  array  for  the  next  bearing  and 
so  on  in  a  pipeline  fashion.  Figure  13  illustrates  the  relationship 
between  the  processes. 

3.  Semaphores 

Communication  between  the  processes  was  accomplished  with 
binary  semaphores.  The  Indata  process  set  the  semaphores  A_data  and 
Bjdata  whenever  it  was  utilizing  either  the  A_indata  or  B_indata  array. 
Similarly,  each  of  the  parallel  Sum_amplitude  processes  set  a  semaphore 
depending  on  the  data  array  it  was  processing.  Thus  for  three  parallel 
Sum_amplitude  processes  six  semaphores  were  utilized:  A_ampl ,  A_amp2, 
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A_amp3,  and  B_ampl ,  B_amp2,  B_amp3.  The  Indata  process  waited  for  each 
of  the  Sum_amplitude  semaphores  for  one  of  the  data  buffer  arrays  to  be 
unlocked  before  it  could  proceed  to  fill  the  buffer  with  new  sampled 
data  from  the  hydrophone  outputs.  Similarly,  each  of  the  Sum_ampli- 
tude  processes  waited  on  the  A_data  or  B_data  semaphore  before  It 
began  processing  the  data.  In  this  way  completely  independent  but 
synchronized  operation  of  all  the  processes  was  achieved.  The  opera- 
tional form  of  the  algorithm  in  Appendix  C  illustrates  the  parti- 
tioning of  the  processes  and  the  use  of  the  semaphores. 

E.   EXPERIMENTAL  DEMONSTRATION 

1 .  Equipment 

To  test  and  demonstrate  the  feasibility  of  near  real-time 
execution  of  the  redesigned  algorithm  an  Intel  lee  800  Microcomputer 
Development  System  (MDS)  [7]  ,  shown  in  Figures  14  and  15,  with  the 
CP/M  monitor  control  program  [s]  was  used.  The  MDS  is  capable  of 
being  used  as  a  partial  system  simulator  to  check  out  software  execu- 
ting from  RAM  in  the  Intellec  MDS.  A  total  of  64  kilobytes  of  memory 
was  available  and  four  Intel  SBC  80/20  Single  Board  Computers  [5]  were 
inserted  into  the  MDS  motherboard  as  shov/n  in  Figures  16  and  17. 

2.  The  Processes 

The  contents  of  the  Trace  matrix,  the  pre-calculated  traces, 
were  computed  using  the  program  of  Appendix  B,  and  stored  in  a  floppy 
disk  data  file  for  subsequent  loading  into  memory.  The  redesigned 
"operational"  algorithm  (Appendix  C)  was  written  as  if  it  were  to  be 
executed  on  a  single  processor,  and  consisted  essentially  of  a  Setup 
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procedure  that  calculated  the  contents  of  the  address  lookup  arrays,  the 
Indata  procedure,  and  the  Sum_ampl1tude  procedure.  This  was  compiled  and 
linked  so  the  program  started  at  location  3000  Hex.  The  choice  of  3000  H 
as  the  initial  address  for  the  program  was  based  on  the  fact  that  addres- 
ses 3000  H  to  3EFF  H  were  available  as  private  RAM  on  each  of  the  SBCs. 
All  the  data  structures  and  variables  such  as  the  semaphores  that  were 
to  be  global  to  all  the  processes  were  located  in  common  memory  above 
4000  H  that  was  accessible  by  all  the  SBCs.  An  empty  dumny  array  of  the 
necessary  size  was  declared  to  ensure  the  assignment  of  common  memory 
addresses  to  the  global  variables. 

Using  the  MDS  and -the  Dynamic  Debugging  Tool  (DDT)  program  [9] 
the  program  was  altered  at  the  machine  code  level  (hex  representation) 
by  removing  the  appropriate  call  statements  and  creating  two  versions: 
one  which  executed  only  the  Indata  procedure,  and  the  other  the  Sum_am- 
plitude  procedure.  In  this  way  the  two  processes  stored  in  the  private 
memory  of  separate  SBCs  could  operate  independently  but  in  a  syn- 
chronized manner  using  the  semaphores  in  common  memory.  The  data  that 
was  operated  on,  as  well  as  all  the  data  structures,  were  available  in 
common  memory. 

The  Sum_amplitude  version  was  then  amended  using  DDT  so  that  all 
the  local  variables  in  the  procedure,  which  the  compiler  had  assigned  to 
memory  locations  in  the  common  area,  were  readdressed  so  they  were 
within  an  SBC's  private  memory.  The  locations  assigned  to  the  dummy 
array  were  utilized  for  this  purpose.  This  then  provided  a  template  for 
three  subsequent  copies  of  the  Sum_amplitude  process  each  of  which  was 
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individually  altered  to  process  a  segment  of  the  range  indices,  and  then 
stored  on  a  disk  file. 

Similarly,  the  addresses  of  variables  that  were  local  to  the 
Indata  procedure  were  amended  so  they  were  within  an  SBC's  private 
memory.  Having  loaded  the  Trace  matrix  from  the  floppy  disk  data  file 
into  the  appropriate  memory  locations  that  had  been  assigned  by  the 
compiler,  the  entire  Indata  version  (which  included  the  data  structures) 
was  saved  on  a  disk  file.  A  subsequent  reloading  of  this  file  into 
memory  caused  the  Trace  matrix  to  also  be  available  and  so  a  separate 
loading  of  the  Trace  matrix  data  file  was  not  necessary.  For  the 
experimental  demonstration  of  the  feasibility  of  the  multiple  micropro- 
cessor concept  the  Indata  process  controlled  the  bearing  parameter  and 
caused  the  entire  algorithm  to  scan  ten  complete  sectors  and  then  stop. 
An  operational  version  would  continue  indefinitely.  As  the  double  buf- 
fering design  resulted  in  two  bearing  increments  being  processed  on  each 
pass  through  the  procedure,  it  was  necessary  to  ensure  that  each  sector 
consisted  of  an  even  number  of  bearing  increments,  i.e.  thirty  vice 
thirty-one.  No  attempt  was  made  to  include  pulse  propagation  time  or 
data  acquisition  time.  The  procedure  that  would  fulfill  these  functions 
was  left  as  a  "stub"  since  it  was  not  required  for  the  demonstration  of 
the  synchronization  between  the  processes. 
3.  Results 

The  capabilities  of  the  MDS  and  DDT  were  use3  to  subsequently 
load  the  Indata  process  into  the  private  memory  of  one  of  the  SBC's,  and 
the  Sum_amplitude  processes  into  the  private  memory  of  each  of  the  three 
remaining  SBCs.   Figure  18  provides  a  block  diagram  of  the  experimental 
configuration  where  process  D  is  the  Indata  process  and  processes  S-j  to 
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Figure  18.   Block  diagram  of  experimental 
configuration. 
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$2   are  the  Sum_ampl1tude  processes.  With  the  processes  executing  simul- 
taneously the  time  required  for  the  amplitude  summation  function  to  be 
performed  on  ten  sector  scans  was  recorded. 

Executing  the  algorithm  with  a  different  number  of  Sum_ampli- 
tude  processes  in  parallel,  while  ensuring  that  all  range  indices  were 
processed,  permitted  an  evaluation  of  the  effect  of  increasing  the  number 
of  parallel  processors.  The  resulting  times  for  one  sector  scan  and 
consequently  the  performance  in  terms  of  scans  per  second  are  given  in 
Table  2  and  Figure  19. 


Processes 

Time  for  one 
(seconds  ±0. 

scan 
.01) 

Scans  per 
Second 

D.  S^ 

6.45 

0.155 

D,  S^.  S2 

3.23 

0.310 

D,  S^,  S2.  S3 

2.15 

0.465 

TABLE  2.  Results  of  the  experimental  demonstration  of  the 
feasibility  of  implementing  the  signal  processing 
algorithm  on  multiple  microprocessors. 


The  correct  operation  of  the  algorithm  together  with  the  per- 
formance times  indicated  that  the  design  was  a  feasible  approach  to 
obtaining  near  real-time  execution  of  the  proposed  signal  processing 
algorithm.  The  fact  that  the  processes  were  stored  in  the  SBCs  private 
memory  and  required  access  to  the  bus  only  when  the  common  data  struc- 
tures and  semaphores  were  referenced  resulted  in  a  minimum  of  bus  usage 
An  examination  of  the  code  for  the  innermost  loop  of  the  Sum_amplitude 
process  (the  summation  of  the  Indata  array  elements  identified  by  the 
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trace)  indicated  that  the  bus  usage  was  7.9^  of  the  loop's  execution 
time.  The  number  of  parallel  SBCs  executing  simultaneously  could  theo- 
retically be  increased  to  twelve  before  bus  contention  would  adversely 
influence  the  overall  execution  time.  Table  2  and  Figure  19  depicting 
the  results  of  the  experiment  clearly  indicate  the  linear  relationship 
between  performance,  in  terms  of  scans  per  second,  and  the  increasing 
number  of  parallel  processors  used.  If  the  number  of  SBCs  operating  in 
parallel  were  increased  to  seven,  plus  one  for  the  Indata  process  for 
a  total  of  eight,  then  a  scan  rate  of  one  sector  per  second  would  be 
achieved. 
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V.   CONCLUSIONS 

For  the  problem  of  locating  and  classifying  objects  buried  in 
marine  sediments  an  underwater  acoustic  imaging  system  has  been  pro- 
posed and  the  thrust  of  this  thesis  has  been  the  design  of  the  signal 
processing  algorithm  that  would  extract  the  information  from  scattered 
acoustic  signals  for  subsequent  operator  display.  Using  a  developmen- 
tal form  of  the  algorithm  it  was  determined  that  the  acoustic  imaging 
system  should  be  a  narrow  beam  system  for  optimum  resolution  and 
imaging  performance.  Since  the  processing  time  was  excessive,  the 
algorithm  was  redesigned  for  implementation  on  multiple  processors  for 
execution  efficiency.  It  was  experimentally  demonstrated  that  this 
design  was  a  feasible  approach  to  obtaining  near  real-time  execution  of 
the  algorithm. 

The  maximum  rate  of  sector  scanning  that  is  possible  for  the  pro- 
posed system  is  slightly  in  excess  of  two  sector  scans  per  second. 
This  was  arrived  at  by  considering  the  pulse  propagation  time  and  the 
data  acquisition  time,  a  total  of  14.5  milliseconds  at  each  bearing, 
as  the  minimum  time  between  pulses.  It  has  been  shown  that  the  use  of 
multiple  microprocessors  can  result  in  an  efficient  implementation  of 
the  proposed  signal  processing  algorithm.  The  number  of  parallel  pro- 
cessors required  to  reduce  the  amplitude  summation  process  to  about  15 
milliseconds  (and  thus  achieve  the  maximum  rate  of  two  scans  per 
second)  would  be  dependent  on  the  microprocessor  chosen. 
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VI.   RECOMMENDATIONS 

The  signal  processing  algorithm  and  its  implementation  on  mul- 
tiple microprocessors,  as  discussed  in  this  thesis,  would  form  the  core 
of  a  special  purpose  processor  for  the  acoustic  imaging  system.  Two 
areas  would  require  development  and  detailed  definition  before  the 
entire  special  purpose  processor  could  be  designed:  data  acquisition 
to  provide  inputs  to  the  processing  algorithm,  and  output  from  the 
algorithm  in  the  form  of  a  display. 

It  is  recommended  that  the  Indata  process  outlined  in  the  thesis 
be  expanded  from  the  current  "stub"  status  to  perform  the  functions 
previously  mentioned;  i.e.  control  of  pulse  transmission  timing, 
sampling  the  hydrophone  outputs,  and  storing  the  results  into  the  two 
buffer  arrays.  It  is  anticipated  that  the  approaches  used  for  the 
algorithm  design  would  be  readily  applicable  to  these  functions.  How- 
ever, the  actual  details  would  be  highly  hardware  dependent  which  would 
require  investigation  prior  to  designing  a  parallel /pipeline  approach  to 
this  portion  of  the  specialized  processor. 

Similarly,  the  interface  between  the  multiple  microprocessor 
system  (performing  the  signal  processing)  and  a  display  processor  would 
require  definition  and  would  be  hardware  dependent.  It  is  considered 
likely  that  the  image  information  after  each  sector  scan  could  be 
stored  thus  forming  a  three  dimensional  representation  in  memory  (the 
third  dimension  resulting  from  the  travel  of  the  underwater  vehicle 
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over  a  section  of  sediment).  This  would  provide  an  opportunity  for  a 
two  dimensional  display  in  either  the  vertical  plane  (such  as  range  and 
bearing)  or  in  a  horizontal  plane  at  some  range.  Furthermore,  the  possi- 
bility of  a  three  dimensional  display  could  be  considered.  It  is  evi- 
dent that  there  is  much  promising  work  to  be  accomplished  in  order  to 
produce  the  special  purpose  processor  for  the  acoustic  imaging  system. 
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